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Abstract 

In this paper, we propose a feature-based method for spectrum sensing of OFDM 
signals from sub-Nyquist samples over a single band. We exploit the structure of the 
covariance matrix of OFDM signals to convert an underdetermined set of covariance- 
based equations to an overdetermined one. The statistical properties of sample covari¬ 
ance matrix are analyzed and then based on that an approximate Generalized Likeli¬ 
hood Ratio Test (GLRT) for detection of OFDM signals from sub-Nyquist samples is 
derived. The method is also extended to the frequency-selective channels. 


1 Introduction 


Cognitive Radio (CR) is emerging as a promising technology for improving the efficiency of 


radio spectrum use in wireless communication systems (IMitola and Maguire! 


1999 1. Spec¬ 


trum sensing (SS) is the most vital task in CR defined as identifying spectrum holes by 
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sensing th e rad i o spec trum and utilizing them without causing interference to primary users 


(PUs) ( 

Havkin et al. 

, 2 

009 

). Of special int 

erest 

in this regard is sensing of OFDM signals 

(Chaudhari et ah. 

2009 

: Axcll and Larsson. 

2011 

). OFDM is one of the most effective mul- 


ticarrier techniques for broadband wireless communications which is employed by many of 
the current and emerging wireless technologies. 


On the other hand, due to the limitations of today’s analog-to-digital converter (ADC) 
circuits which cannot support very high bandwidth and need excessive mem ory and p ro- 


hibitive energy costs for implementing digital signal processing systems (jCohen et al 


20111 ), 


it may be very costly and even impractical to sense the signal based on Nyquist-rate sam¬ 
ples. This has motivated researchers 


Tian et al.. 

2012 


Polo et ah. 

2009; 

2011; 

Cohen et al. 

201.1; 

Sun et ah 


(to 

lishali 

and Eldar. 

2011 

Tian and Giannakis. 

2007; 

Tian. 

2011 

Leus and Tian, 

2011; 

Ariananda and Leus. 

2012; 

Rebeiz et ah, 

2012 

). But to the best of our knowl- 


edge, there has been less efforts targeting the detection of OFDM signals from sub-Nyquist 
samples. In this paper, we propose a new feature-based approach for sensing an OFDM 
signal occupying a single band from sub-Nyquist samples. 


Related research and previous work: The problem of OFDM sensing using second- 


order statistics has been alrea 

dv studied in. e.g., ( 

Chaudhari et ah, 

2009; 

Axell and Larsson. 

2011; 

Bokharaiee et ah. 

2011; 

Al-Habashna et ah 

2012 

). All of these methods in some way 


exploit the correlation induced by CP to sense the presence of OFDM signal, but perform the 
detection based on Nyquist rate samples which, as discussed earlier, might need expensive 
ADCs espeially for wideband signals. 


During the past several years, the problem of spectrum sensing from sub-Nyquist samples 
has attracted a lot of attention. While many of the approaches propos ed so f a r re ly on the 


sparsi ty in_t 


(120071 1: 


re frequen cy d omain 


Tianl (120081 1: 


Polo et al 


arising from spectrum underutilization 


Tian and Giannakis 


(12009); 


Sun et al. 


(12013! ). there has also been methods which 
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Leus and Tian 


tall); 


Cohen et al. 


(201l|); 


Rebeiz et al. 


do not necessarily need sparseness in the sr 

>ectrum to work 

Tian et al. 

(2012 

); 

Tian 

Leus and Tian 

(20 

LD 

; Ariananda and Leus 

(2011) 

Razavi e 

al. 

(20lj); 

Cohen et al 

Rebeiz et al. 

( 

>012 

). For instance, the met 

Lods proposed in 

Tiai 

et al. 

(2012 

): 

Tian 


( 2012 ); 


Razavi et al. 


( 2011 ); 


( 2011 ); 


(mil); 


( 20131 ) are 


based on recovering the Spectral Correlation Function (SCF) of the signal from sub-Nyquist 
samples. All of these methods need sparsity in the SCF for de tecting the s ignal. Since the 


SCF matrix of the Nyquist-sampled OFDM signal (see, e.g., (ITian et al. 


2012 


Equation 


(10))) is nonzero everywhere, then it is not possible to use these methods for recovering 
the SCF of OFDM signals from sub-Nyquist samples as they need sparsity in the SCF for 


recovering it using Compressive Sensing (CS) algorithms or converting 


set o f equati ons to an overde 


( 2011 1 and (Tian et al. 


ermined one. The methods proposed in 


he underdetermined 


Ariananda and Leus 


2012 . Section IV) reconstruct the Power Spectral Density (PSD) 


from sub-Nyquist samples and then decide on the presence of sign a l based o n the recov¬ 


ered PSD. Similar to conventional energy detection methods 


Yucek and Arslan (20091 ). the 


main challenges with these meth ods are s electing the threshold and poor performance in 
low signal-to-noise ratios. In 


Razavi et al. 


( 20131) . we detected an OFDM signal over a sin¬ 


gle band from sub-Nyquist samples. Although this method assume that the whole band is 


occupied 

signature 


w the OFDM signal, but we inj e ct th e sparsity to SCF using a cyclostationary 


Adrat et al. 


(2009); 


Sutton et al. 


( 20081 ) embedded in the signal to assist with the 


detection procedure. This migh ; 


work identification 


Razavi et al. 


)e of help in the problem of rendezvous and cognitive net- 


( 20141 ). but it cannot be used in detecting primary users, as 


primary users usually do not care about secondary users and therefore do not assist them 
via injecting a signature in their own signals to facilitate the spectrum sensing task. 


Our contribution: In this paper, we propose a new feature-based approach for sensing 
OFDM signals from sub-Nyquist samples over a single band entirely occupied by the sig¬ 
nal. In other words, we do not assume any sparsity in the spectrum. Instead, the method 
exploits the features of covariance matrix of OFDM signals stemming from the insertion of 


3 











































































































cyclic prefix (CP). These features, as we will show later, will help us to convert the un¬ 
derdetermined set of two-dimensional (2-D) equations which relates the available covariance 
matrix of sub-Nyquist samples to the unavailable covariance matrix of Nyquist samples, to 
an overdetermined one-dimensional (1-D) one. Then, based on the statistical properties of 
the sample covariance, an approximate GLRT-based detector is derived. The method is also 
extended to the case of frequency-selective channels. 

Paper organization: The rest of the paper is organized as follows. In Section [2] the 
system model over Gaussian channels is given and the problem is formulated. In Section [3] 
the relevant statistical properties of sample covariance matrix are studied and a covariance- 
domain linear system is derived. The approximate GLRT-based detector is introduced and 
formulated in Section [4j In Section [5] we extend the results to the frequency-selective channel 
case. Finally, we study the performance of the proposed method by simulation experiments 
in Section |6] Conclusions are drawn in Section |T] and some details of the derivations are 
given in the Appendix. 


Notations and Mathematical Preliminaries: Throughout this paper matrices and vec¬ 
tors are denoted by capital and small boldface letters, respectively. = denotes the equality 
and = denotes the definition. E is reserved for statistical expectation operator, Cov(x, y) 
represents the covariance matrix between random vectors x and y, and <g> denotes the Kro- 
necker product. Ip and 0 pq represent, respectively, PxP identity matrix and PxQ full-zero 
matrix. For an arbitrary M x N matrix A, [A]jj denotes the (i,j)-th entry of the matrix 
and ap % — 1, 2,..., N denotes its i-th column. vec(A) is the vectorization of A defined as 
vec(A) = [a^f,..., a Jf] T . Furthermore, if M = N, then vech(A), known as half-vectorization 


of A, is the Jv ^ +1 - > x 1 vector which is obtained by column-wise stacking of only 


he elements 


2008 


on and below the main diagonal. It can be proven that for an N x N matrix A (jSeberf 
Chapter 11) 

vech(A) = f2jvvec(A), (1) 
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where 


T (0) 

l n 

®N,N 

Oat,at 

Oat,at 


N-1,N 

I (1) 

Oat-1,AT • 

• Oat-1,at 


N-2,N 

Oat-pat 

t(2) 

l n 

• Oat-2,AT 

(2) 

Ol,N 

Oi,at 

Ol,AT 

t (AT-1) 

i-N 



where 1^ denotes an identity matrix whose fi rst i rows are discarded. Remark that I® — In- 


Moreover, if A is a symmetric matrix, then flSeber . 


2008 


, Chapter 11) 


vec(A) = rVvech(A), 


(3) 


where T is an (IV 2 ) x (N(N + 1)/2) matrix with entries ((i — 1 )N + j, (j — 1 )(N — j/2) + i) 
and ((j — 1 )N + i, (j — 1)(A— j/2) + i) for 1 < j < i < N equal to 1 and the rest of elements 
equal to zero. Remark that = Ijv(jv+i)/ 2 - 


Furthermore, for any th ree ar 


following two equalities ( Sebe r 


j itrary matrices A, B, and C of suitable sizes, we have the 


2008 


Chapter 11) 


vec(ABC T ) = (C<g> A)vec(B). 


(4) 


and 

ABC T = ^^[B], 3 a,cJ, (5) 

i j 

where a.j and Cj are i-th and j-th columns of A and C, respectively. 
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2 System Model over Gaussian Channels 


Consider a secondary user (SU) with the goal to detect the presence of an OFDM signal. 
Denoting the received signal by x(t), we can formulate the problem as deciding between the 
following two hypotheses 


{ "Ho : x(t) = nit ) 

( 6 ) 

Hi : x(t) = s(t) + nit) 

where s(t) is an OFDM signal and n(t) ~ CAf(0, cC) . If the number of subcarriers is large 
enough, from central limit theory we have s(t) ~ CA7(0, of). Suppose next that we have 
sampled the signal x(t) at a sub-Nyquist rate to collect compressive samples z{t). Matrix- 
wise, this can be described a|] 


z (k) = Ax(fe), k = 1 , 2,..., Nb, 


(7) 


where z(k) = [z[fcM], z[kM + 1],..., z[kM + M — 1]] T consists of sub-Nyquist (compressive) 
samples, x(fc) = [x[kN],x[kN- 1-1],..., x[kN+N — 1]] T consists of (unavailable) Nyquist sam¬ 
ples, Nb is the total number of taken frames, and A is the M x N real-valued measurement 
matrix, p = < 1 is called compression (or downsampling) ratio. 


Now, the task of a sub-Nyquist OFDM detector is to decide whether the OFDM signal exists 
or not, based on sub-Nyquist measurements {z(fc)}^f r 


It can be easily seen that R 2 = AR a .A 7 where the N x N matrix R s = E(x(/c)x^(/c)) 
and M x M matrix R 2 = E(z(fc)z i/ (fc)) are covariance matrices of x and z, respectively. Let 
us denote the useful symbol length of the considered OFDM signal by T, the cyclic prefix 


x We remark here that there are two main strategies for sub-Nyquist sampling of a signal iMishali ct~al 


(2011): Random Demodula tor (RD), also known as Analog-to-Informatio n Converter (AIC) . Tropp et al 
( 2010 ): Kirolos et al. ( 200fil) . and Modulated Wideband Conve rter (MWC) Mish ali and Eldail (|2010h . Both 
of these strategies can be mathematically formulated as in ([7]) I Mishali et al.1 (|201lh . 
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length by T cp and the total OFDM symbol length by T sym = T + T cp which are all assumed 
to be known to the cognitive user. Assume that the basic sample duration is normalized to 
1. Then by setting N = T sym , matrix R x of the noisy OFDM signal (i.e. x(t) under "Hi) can 
be written as 


a 2 s + a 2 if i = j, 




err 


if \i-j\ = T, 
otherwise. 


( 8 ) 


On the other hand, under "H 0 we have 




if * = j, 

0 otherwise. 


(9) 


To encompass the descriptions of R x under TLq and "Hi, we present it by 


Rx = tqIn + T s A, 


( 10 ) 


where In is the N x N identity matrix and A is defined as 


[A]y ^ 


A , 1 H\i-j\=T, 


0 otherwise, 


( 11 ) 


and then distinguish between T-L 0 and 1~L\ as 


T-io ■ (to,t s ) = {al, 0 ), 

■ (t 0 ,Ts) = {al + a 2 s ,a 2 s ). 


( 12 ) 
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Since in both cases of (fT2|) we have r 0 = t s + therefore (TT2|) can be simplihed as 


j H 0 : t s = 0 , 
| Hi : r s 0 . 


( 13 ) 


In other words, the problem of detection of OFDM signal can be re-expressed as identifying 
whether in the general description of R x in (ITU the parameter t s is zero or not. 

In practice, covariance matrices are not readily available and hence we substitute them by 
sample covariance matrices R x = Ylk=i *-(k)x- H (k) and R z = ^ z(k)z H (k), where 
Nb is the number of frames used for the computations. It is easy to see that 


R z = AR X A 7 . 


(14) 


In the next sections we will show how the sub-Nyquist sample covariance matrix R z can be 
used to perform the hypothesis testing in (J6J) based on the observation made in (1I3|) . 


3 Covariance-Based Linear Equations for OFDM De¬ 
tection 


The sample covariance matrix R x can in general be written as the sum of covariance matrix 
R x and an error term W stemming from the finite-sample effects. From now on, we call W 
the finite-sample noise and express this as 


R r — Rx + W. 


(15) 


It is easy to verify that E X (R X ) = R x under both Ho and Hi , and therefore conclude that 


E X (W) = 0. Furthermore, from (IGoldbergerl. 


1991 


, Chapters 9-10) we know that the entries 





of W have asymptotic Normal distribution. The following th eorem t hen states th e cov ariance 


of entries of W. Parts of this theorem can also be found in 


Axel l and Larssonl f2 01 il l. 


Theorem 1. Suppose that Wij and w P)q are two arbitrary entries of W below or on main 
diagonal (i.e. i > j, p> q). Then 


if if, j) = (p,q) 


®(wijw Ptq \n 0 ) = \ 2Nb 

0 otherwise , 


(16) 


where r 0 = of, and 


T o 

Ng 


E(w itj w Piq \Til) = < 


if(i,j) = ( P,q),i~j = 0, 

T?+r " */ (i,j) — ( P, q),i~J=T, 

if (hj) = ( P,q),i~j i {0,T} 


r 2 +r 2 

~s +~Q 

2Ng 


T o 

2N b 


tfi=j=P=q+T 


otherwise , 


(17) 


where To = + of and t s = a 


A 


Proof. Proof is deferred to Appendix. 


□ 


Inserting next (1T5T) in (IT4l) yields 


R z = AR t A J + AWA J 


(18) 


Notice that since R x and A are both real-valued, the signal part of (I18p is real-valued and 
the imaginary part only includes the effect of finite-sample error W. In fact, we can simply 
change (TT8li to a real-valued equation by keeping only the real part of R = and throwing away 
the uninformative imaginary part. To avoid introducing extra notations, from now on we 
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assume that (1T8|) represents a real-valued matrix equation. 

The problem of sensing OFDM signal can then be restated as testing whether in 2-D linear 
model (pgp we have R a . = To I or R a = rol + r s A for some nonzero unknown parameters To 
and t s . The first problem in this regard is that (1T8T) represents an underdetermined linear 
system of equations. To solve this issue, we first apply the vech operator to (HHD and use 
o. a, a and PD to obtain 

r z = n M ( A ® A)(r 0 vec(I) + T s vec(A)) + Cl M ( A ® A)Rvw, (19) 

where r~ = vech(R 2 ) and w = vech(W). From (J4J) and (J5l) it is easy to verify that (A <g> 

N N 

A)vec(I) = ® a *) an d (A <8> A)vec(A) = (a<8) a^ + a* (g) aj_y) and therefore 

i =1 i=T+l 

rewrite (HD as 


r z = T s b s + T 0 b 0 + v 

= B0 + v (20) 

N N 

where b 0 = fl M Z)( a * ® a *) and b s — I] (a,_ T ® a* -fa,® a *_ T ), B = [b 0 b s ], 

i= 1 i=T +1 

0 = [to t s ] t , v = Kw is the additive noise vector, and K = CIm {A <g) A) I'm- 

Now, the problem of OFDM signal detection can be expressed as testing whether in the 
classical linear model ()20l) the second entry of parameter vector 9 is zero or not. We formulate 
this in the next section. 

We remark here that since in linear regression equation (1201) the number of unknown param¬ 
eters is two and the number of equations is M(M + l)/2, any M > 2 is theoretically enough 
for estimating r s and performing the test. 
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4 Approximate Generalized Likelihood Ratio Test 


To derive the GLRT-based detector, we first notice that the noise vector v in (l20]i is not 
white. Therefore the first step is to whiten the noise by multiplying both sides of (1201) by 
£“ 1/2 , where 

= K£ w K t , (21) 

and is the covariance matrix of w = vech(W). Remark that as proved in Theoremdl T, w , 
and therefore £„, does not expose the same expression under "Ho and "Hi. Hence, to simplify 
the test, in this section we formulate an GLRT-based detector based on approximating 
£ t , with a hypothesis-independent covariance matrix which compromises the properties of 
covariance matrices under the two hypotheses. 

We first adopt a hypothesis-independent approximation for the entries of T, w as follows 


E (WijWp^g) 


2c if (i,j) = (p,q),i~j = 0, 
< c if (i,j) = (p,q),i~j ± 0 
0 otherwise, 


( 22 ) 


where c is an unknown constant. Remark that W,(wijW Piq ) in ([22]) shares properties with 
both E(w i jW Ptq \'H 0 ) in (IHUl and hi (lT7[k like (fTRll it is zero for (i,j) ^ (p,q), 

but then for (i,j) = ( p,q ) it is like (ITT]) with t s r 0 . Equation (122|) then implies that S w 
has the following form 

S w = C A, (23) 

where A is an N( - N + 1) diagonal matrix with diagonal elements in locations 1, 1 + (N), 
1 + (N) + (N — 1), 1 + (N) + (N — 1) + (N — 2), ... being equal to 2 and the rest are equal 
to 1. In fact, the diagonal elements having value 2 correspond to the first case of (1221) and 
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the rest correspond to its second case. The third case of (122]) implies that H w is diagonal. 
Now, multiplying both sides of (j20[) by T = (KAK T ) -1 / 2 yields: 


r 2 = B0 + v, 


( 24 ) 


where r z = Tr z , B = TB, and v = Tv. Based on this, the problem of OFDM detection can 


be expressed as doing hypothesis testing problem (1T3l) for linear system (j24l) in which the 
additive noise has distribution v ~ Af(0, cl) with some unkn own variance c. The GLRT- 


based detector for this problem can then be written as flKavl . 


19981 . Theorem 9.1) 


T(f z ) 


A 


(N' 


2)(D0 1 ) t [D(B t B)- 1 D t ]- 1 (D0 1 ) «i , 

- \ 7 

r^I — B(B T B)- 1 B r jr 2 «o 


(25) 


where N' = j ' / 0^+ 1 ) denotes the number of distinct equations, D = [0 1], and 0 1 = 
(B T B) _1 B T r z is the maximum likelihood estimate of 6 under "Hi- Furthermore, the proba¬ 
bility of false alarm (Pfa) and probability of detection (Pd) of the approximate GLRT-based 
detector are given by 


Pfa - QF 1N ,_ 2 (y), 

(26) 

Pd = Qf 1jv ,_ 2 (a)(7 , ) ; 

(27) 


where the noncentrality parameter takes the form 


(D6>i) t [D(B t B)- 1 D t ]- 1 (D0 1 ) 

c 


(28) 


where 6\ is the true value of 6 under Hi. 
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5 Extension to Frequency-Selective Fading Channels 


In this section we extend the results obtained in the previous sections to the case of wideband 
frequency-selective channel. Assume that the multipath channel between the PU transmitter 
and the SU receiver has the model H(z) = h 0 + h\Z + ... + h L z~ L . Putting the tap weights 
in (L + 1) x 1 vector h = [ho, hi,, h k ] T and denoting the k -th OFDM block of length 
T S ym = N as s fc = [sfc(l), Sfc(2),. .. ,s k (N)] where s k {n) = s[(k - 1 )N + n], we have 


x fc = S k h + n, k = 1,2,..., N B , 


where S k is a Toeplitz N x [L + 1) matrix defined as 


Sfc(l) 

Sfc-i(iV) 

Sfc-i(Af - 1) 

Sfc(2) 

Sfc(l) 


Sfc(3) 

S k ( 2) 

Sfe(l) 


s k ~i(N — L + 1) 
s k- i(AI — L + 2) 
s k ~i(N — L + 3) 


Sfc(A^) s fc (V-l) s k (N-2) 


s k (N-L ) 


Denoting the i-th column of Sfc by we can write 


"P J u lj 


(29) 


(30) 


(31) 


where J d and J u are respectively down-shift matrix and up-shift matrix whose (i,j)-th 
elements are [Jd]i,j — $i,j+i and [J u ]ij = fii+i.j where 8 it j denotes the Kronecker delta. 
Remark that J d = J 
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From (ISeberl . 


2008 


, Theorem 21.6) and since E(Sfc) 


Oat,l+i we can next write 


Rx = E(S fc hh ff Sf) + all 

L L 

= EE hih* Cov(s M , s fc j) + ctqI. (32) 

i =0 j =0 


From (l3Tlh Cov(sfc i j, Sfej) can be computed as 


Cov(s fcii , 8*,^) 


= ^ (j;<a+1) +1) 


(33) 


Denoting c, ;j = vec^J^(A + I)J{ + *(A + I)J^ , we can write 


vec(Rx) 


L L 

h i h *3 C hi + CT 0 vec ( I ) 

4 = 0 j =0 


[C s , vec(I)] 


Un 


(34) 


where C s is dehned as the N 2 x (L + l) 2 matrix whose (j(L + 1) + i + l)-th column is c h j 
(remark that 0 < i, j < L) and h is dehned as the vector whose (j(L + 1) + i + l)-th entry 
is a 2 hih*. 

Putting (jHD hi (TT8T) and vectorizing the resulting equation yields 


r z = B m e m + v (35) 

where B m = f2 M (A <g) A)[C S ,vec(I)], = [h 7 ,ao] T , and v = Q M (A <g) A)Fjvvech(W). 
Then the problem of OFDM detection can be expressed as performing the following test for 
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the linear model of (135|) 


j Ho : D m,0m = 0, 

^ H\ : D m G m 7^ 0, 


( 36 ) 


where D m = [I( i+ i) 2 , 0 (i + i) 2 jl ], The rest of the test formulation is similar to the Gaussian 
case discussed in 01 


6 Simulation Examples 

In this section we study the performance of the proposed sub-Nyquist OFDM detector by 
simulation examples. For simplicity, we consider an OFDM system with IFFT size 32 which 
means T = 32. The cyclic prefix length is chosen as T cp — Tj 4 = 8 and subcarrier symbols 
are assumed to be drawn from a 16-QAM constellation with unit energy. The elements of the 
measurement matrix A are drawn from A/"(0,1) and its columns are normalized to have unit 
norm. We emphasize that these chosen values are just examples for carrying out numerical 
simulations, and are not as such related to the fundamentals of the derived detectors in any 
way. 

Before proceeding with the simulation examples, we briefly discuss the parameters which 
affect the performance based on the formulation of the problem given in previous sections. 
We first recall that the linear equation used for detecting the OFDM signal is (fl8jl . Therefore 
the parameters which affect the detection are those affecting the solution of this 2-D linear 
regression problem. The first parameter is the number of independent equations in (fTSlh i.e., 
M(M + l)/2. For a fixed Nyquist frame size A, this parameter is uniquely specified by the 
compression ratio p. The second and third parameters are those which affect the power of 
finite sample error W. Formulas clearly show that, for a given signal power o 2 , this 

is determined by iVg; the number of blocks taken for computing R,, and a 2 ; the variance 
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of noise, or equivalently, signal to noise ratio. The first two experiments in this section are 
devised based on the above observations. Besides, the exactness of the approximation we 
adopted in Section [4] will be studied by simulation results. As mentioned in Section □] there 
is no method in the literature specifically designed for single-band OFDM detection from 
sub-Nyquist samples. Therefore , to co mpare our method with some existing methods, we 


choose the method proposed in (ITian et al. 


2012. Section IV) which recovers the PSD of a 


general stationary signal from sub-Nyquist samples. 


The third experiment studies the performance of the method over frequency-selective chan¬ 
nels. 


6.1 Influence of the compression ratio 


In the first experiment we study the effect of compression ratio on detection performance. 
For this experiment, the detector deploys Ng = 100 OFDM symbols. Figured] illustrates the 
probability of detection of the proposed method as a function of signal-to-noise ratio (SNR) 
for different compression ratios p G {0.2, 0.4, 0.6, 0.8,1} when the threshold is set to obtain 
Pfa = 0.05. To check also the exactness of the approximation adopted in Section [4] we 
compute the threshold y' both from the true Pfa calculated from the simulated data in the 
absence of OFDM signal (solid lines with circle markers) and from formula (1261) ( dotted lines 
with diamond marker). The results of the PSD recovery method introduced in Ti an et a l. 


(120121) are shown by dashed lines with square markers. For this method the threshold has 


been computed from the true Pfa- 


The first observation from [T] is that, as we expect, the performance of the proposed method 
enhances with increasing the compression ratio. Interestingly, with heavy compression ratio 
of 0.2, the detection probability is still above 95% at an SNR of 0 dB. Notice that since the 
number of independent covariance-based equations is of order 0(M 2 ), a compression ratio 
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of p = M/N = 0.2 implies that the effective compression ratio (in covariance domain) is 
only about 0.2 2 = 0.04 or 4%. Also, as it can be observed, the proposed method provi des a 


much better performance compared to the general PSD recovery method of ( Tian et al. 


2012 


Section IV). This can be seen also from Figure [2] where the Receiver Operating Characteristic 
(ROC) curves of both methods have been depicted. 


Furthermore, the unnoticeable difference between the curves with threshold computed from 
true Pfa and the ones with threshold computed from (j26j) in Figure [1] verifies the exactness 
of the approximation adopted in Section [4] 


6.2 Influence of the number of blocks 


In the second experiment we study the effect of number of blocks taken for computing the 
covariance matrix, Nb, on detection performance. As it can be seen from Figure [3] increasing 
Ns improves the performance of the detector. In fact, in asymptotic case when Nb —> oo 
we have R x = DC Si=u( x (0 x (0’ ff ) = E( xx ^) = Re which means W = 0. Similar to the 
previous experiment, it can be also observed that the performance of th e proposed method 


shows significance improvement over the PSD recovery method of 


Tian et al. 


(20121). This 


can be seen also from Fiure[4]where the ROC curves of both methods have been illustrated. 


Again, the unnoticeable difference between the solid lines (true Pfa) and their corresponding 
dotted lines (target Pfa) in Figure [3] confirms the exactness of the approximation we used 
in Section [ij 


17 








6.3 Performance of the proposed method over frequency-selective 
channels 

The third experiment studies the effect of frequency-selective channels on the performance of 
the proposed method. The method is applied on three different channel lengths with identical 
energy; i.e. o I hi 1 2 is identical for all L = 0,1,2. The OFDM signal characteristics as 
well as the probability of false-alarm are set to the same values as in the two previous 
examples. As it can be observed from Figure [5] the frequency selectivity slightly deteriorates 
the performance but reliable sensing can still be achieved even at low SNRs. 


7 Concluding Remarks 

A method for sensing OFDM signals from sub-Nyquist samples was proposed. The proposed 
method exploits the unique characteristics of the covariance matrix of OFDM signal to 
perform the spectrum sensing task. Based on the statistical properties of the covariance 
matrix we developed an approximate GLRT-based detector. The proposed method was 
also extended to the case of frequency-selective channels. The simulation results verify the 
theoretical observations. The results also illustrate that highly efficient sensing can still be 
obtained, in terms of probability of detection and false alarm, despite of low compression 
ratios and low SNRs. This can open up new possibilities for sensing devices with low-cost 
analog hardware and A/D interface deploying sub-Nyquist observations. 
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Appendix: Proof of Theorem U 


The proof of is trivial. Here we only prove (1171) . Let us define r\_j = Real^[R x ]jjJ, and 
denote the real part of Xi by x t and its imaginary part by x t . Since 


hi 


-i 7V B 

B 1=1 

, N b 

B l _i 


we can write 


Nb N b / 

n,jr P , q = A?E E ([^(0^(0 + 

i=i i'= i' 

[x p (l')x q (l') + ip(/ , )xg(/')]) 

= t x + t 2 , 


(37) 


where T\ and T 2 are defined as follows 

Ng Ng , 

T i - iv^E E + Xi(i)xj(i)] 

i=ir=i v 
i'^i 

[x p (l')x q {l') + x p (l')x q (l')]^, 
[x p (l)x q (l) + x p (l)x q (l)}^ 


The reason for partitioning (1371) to Tj and T 2 is that in T\ the two multiplicative terms 
inside the double summation are independent since they belong to different frames l and l', 
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while in T 2 they belong to the same frame and depending on indices i, j, p, and q might be 
dependent. Due to independence of multiplicative terms, E(T 1 |'H 1 ) is easy to compute as 


Nb-1 ^2 
N b 0 


e(Ti| ni) = 


Ng -1 
N b 


r 0 r s 


N b ~ 1 ^2 
N b ‘s 


0 


if ( i-j,p-q ) = (0,0) 
if (i-j,p-q)e{{0,T),(T,0)} 
if (i-j,P~q) = (T,T) 
otherwise. 


(38) 


Furthermore, after some more scrutinized manipulations, E(T 2 |?^i) can be expressed as 




E(T 2 |'H 1 ) = 




N b T 0 

if 

ihj) = 


O 

1 

2fe( r o+ 3 D 2 ) 

if 

ihj) = 

(p,q),i 

-j = r, 

2ife r 0 

if 

(hj) = 


-j^{o,t}, 

j4 r ° Ts 

if 

i=j = 

P = g 4- 

T. 

-Dr 2 

n b t o 

if 


-q) = 

(0,0 

-L r 2 

iVs D 

if 

(* — 3iP 

-q) = 

(T,T),i^p 


if 

(* — 3iP 

-q)e 

{(0,T),(T,0)}, 




l 

^P 

0 

otherwise. 




(39) 
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Summing (138|) and (l39|) will yield then 




= 


r 2 i A. 

T o + n b 

if (i,j) = ip, q)c — j = o, 

2 , T 0+ r s 

' S " r 2 Ng 

if (i,j) = (p,q),i~j = T, 

T 0 

2Ng 

if (i,j) = i {0 ,T}, 

Ws + ^ 

if i = j= p = q + T, 

T 0 

if (j- fp-j) = (0,0), i yip 

A 2 

if (i-j,P~q) = (' T,T),iyIp 

r 0 r s 

if (i~j,P~q) e {(0 ,T), 


,(T,0)},iyIp 

0 

otherwise. 

) from (140|) and remarking that Wij and r hJ 


(40) 


variances, (Tl7|) is concluded. 
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SNR (dB) 


Figure 1: Probability of Detection versus SNR for various compression ratios p. The number 
of blocks is Nb = 100 and the probability of false-alarm is PpA = 0.05. The solid curves with 
circle markers correspond to the case where the threshold has been determined based on 
the true Pfa computed from the simulated data. The dotted curves with diamond markers 
correspond to the case where the threshold has been determined from formula (l26lh The 
dash ed curv es wi th square markers correspond to the PSD recovery method introduced in 
Tian et al.l (20.12). Color conventions are as follows. Blue: p = 0.2, red: p = 0.4, black: 
p = 0.6, magenta: p = 0.8, green: p = 1. As i t can be seen the proposed algorithm 
clearly outperforms the method introduced in Tian et al. ( 20121 1. Furthermore, the very 
small distance between the solid curves with circle markers and the dotted curves with 
diamond markers shows the exactness of the approximation we adopted. 
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Figure 2: Probability of Detection versus SNR for various compression ratios p. The number 
of blocks is Nb = 100 and the probability of false-alarm is PpA = 0.05. The solid curves 
correspond to the case where the threshold has been determined based on the true Pfa 
computed from the simul ated da ta. The dashed curves correspond to the PSD recovery 
method introduced in Ti an et al.l (2 012 ). The color conventions are as in Figure [0 As i t can 


be see n the proposed algorithm clearly outperforms the method introduced in Tian et al. 

(l2012h . 
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SNR (dB) 


Figure 3: Probability of Detection versus SNR for various number of blocks Nb- The 
compression ratio is p — 0.4 and the probability of false-alarm is Pfa = 0.05. The solid 
curves with circle markers correspond to the case where the threshold has been determined 
based on the true Pfa computed from the simulated data. The dotted curves with diamond 
markers correspond to the case where the threshold has been determined from formula 
(I20D . The dashe d cu rves with square markers correspond to the PSD recovery method 


introduced in Ti an et al 
Nb = 400, black: 


(2012). Color conventions are as follows. Blue: Nb = 100, red: 
Nb = 700, magenta: Nb = 1000. As it can be see n the proposed 

” f|2012h . 


algorithm clearly outperforms the method introduced in Tian et al . (2012). Furthermore, 
the very small distance between the solid curves with circle markers and the dotted curves 
with diamond markers shows the exactness of the approximation we adopted. 











Figure 4: Probability of Detection versus SNR for various number of blocks Nr. The 
compression ratio is p = 0.4 and the probability of false-alarm is PpA = 0.05. The solid 
curves correspond to the case where the threshold has been determined based on the true 
PfA computed from the simu late d data . The dashed curves correspond to the PSD recovery 
method introduced in iTian et al.1 ( 2012 ). The color conventions are as in Figure [3l As i t can 


be see n the proposed algorithm clearly outperforms the method introduced in Tiaii et ah 

(l2012h . 
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Probability of Detection 



Figure 5: Effect of frequency-selectivity of the channels on the performance of the proposed 
method. The channel length is L + 1 and p is the compression ratio. 
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